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In this paper we propose a framework allowing to compute the effective interactions between 
two anisotropic macromolecules, thereby generalizing the DLVO theory to non spherical finite size 
colloids. We show in particular that the effective interaction potential remains anisotropic at all 
distances and provide an expression for the anisotropy factor. We then apply this framework to 
the case of finite rod-like polyelectrolytes. The calculation of the interaction energy requires the 
numerical computation of the surface charge profiles, which result here from a constant surface 
potential on the rod-like colloids. However, a simplified analytical description is proposed, leading 
to an excellent agreement with the full numerical solution. Conclusions on the phase properties of 
rod-like colloids are proposed in this context. 



I. INTRODUCTION 

The DLVO theory, named after Derjaguin, Landau, 
Verwey and Overbeek [1], is one of the most influential 
and still very important description of charged colloidal 
suspensions. It has been developped more than fifty years 
ago to rationalize the stability of lyophobic colloidal sus- 
pensions. One specific prediction of the DLVO theory is 
the far-field pair potential between two spherical colloids 
of like radii a which, within a linearization approxima- 
tion, takes a Yukawa form : 

Z 2 e 2 f exp[K D a]\ 2 exp(-K D r) 

U 12 {r) = — , (1) 

Aire \ 1 + kd<i ) r 

where Z is the valence of the object, e the elemen- 
tary charge and up denotes the inverse Debye screen- 
ing length. The latter is defined in terms of the micro- 
ions bulk densities {p a } (with valencies {z a }) 
Att£b J2 a Pa z a- At the level of a dielectric continuum ap- 
proximation for the solvent with permittivity e, the Bjer- 
rum length Ib is defined as Ib = e 2 /(AirelcsT), where 
ksT the thermal energy: Ib = 7 A for water at room 
temperature. Note that the Debye screening factor, kd, 
does characterize the decay rate of the interaction poten- 
tial in the far field region, providing therefore an exper- 
imental measurement of the screening factor from inter- 
action force measurements (see eg. [2]). 

However, in the colloid world, the spherical shape is 
not the rule and many macromolecules are intrinsically 
very anistropic : rod-like or ribbon-like shapes (DNA 
molecules, TMV or fd virus, V2O5 ribbons, Boehmite 
rods, etc.) [3-7], disk-like shapes (eg for clays, as laponite, 
bentonite, etc.) [8-11]. Since the seminal work of Lang- 
muir on bentonite clay particles published in 1938 [12], 



these systems have been the object of considerable atten- 
tion, in particular in the context of orientational phase 
transitions (such as isotropic to nematic I-N, etc.) [15]. 
From the theoretical side, these transitions were first 
adressed by Onsager [13], who showed that the nematic 
phase was stabilized at high density by purely entropic 
effects. The extension to charged rods has been reconsid- 
ered more recently by Stroobants et al. [14], showing that 
the electrostatic interaction between the polyelectrolytes 
lead to a twisting effect which enhances the concentration 
at the I-N transition. The picture of Onsager reproduces 
correctly the experimental results for highly disymmet- 
ric particles, such as TMV or fd viruses [3,4]. However in 
many anisotropic systems, a gelation occurs before any 
I-N transition [6,8,9]. According to the DLVO theory, 
gelation is usually assumed to result from the presence 
of van der Waals attraction between the macromolecules, 
which overcome at high salinity the double layer repul- 
sion. However, the origin of gelation in many rod- and 
platelet- like systems remains quite obscure [6,9,16]. The 
" gel" denomination is also misleading in some cases since 
the texture of the "gelled" system may be closer to a 
glassy like phase, in which the orientational and transla- 
tional degrees of freedom are frozen [16,17]. The origin 
of such a glass-like transition is still under debate. 

In this paper, we shall stay at a more "microscopic" 
level and consider the effects of anisotropy on the inter- 
action between two macromolecules, much in the spirit 
of the DLVO approach. One specific question we raise is 
the following. We consider two anisotropic particles, sep- 
arated by a "large" distance (i.e. a distance r larger than 
their typical dimension a). Can the electrostatic inter- 
action between these two individual objects be modelled 
by the previous DLVO result, i.e. is the anisotropy lost 
for large distances ? This is of course the case in the 
absence of salt [20]. Does this result generalize with an 
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electrolyte ? 

Before delving into the details, let us first consider 
a much simpler problem, namely that of two identical 
charges q, with positions z = ±a/2 along the z axis (a 
fixed) and embeded in an electrolyte (1): what is the 
electrostatic potential <j) created by these two charges at 
large distances ? 

Naively one would expect that the anisotropy is lost 
for large distances (i.e. distances larger than the size 
a of the object, or larger than the Debye length 1/k.d) 
and the potential should reduce to its Yukawa form 
$(r) = .-|L ex P (-tt D r) _ But thig ig actually not t he case 

! This can be understood by computing -within a lin- 
ear Debye-Huckel like theory- the potential at large dis- 
tances in the x and z directions : along the axis x, one 
gets as expected 4>(r) ~ ^ cx P(~ K P r ) to lowest order 
in a/r; but on the z axis, one gets at the same order 

$(r) _ 2 g cosh(^a) cxp(-K D r) _ Therc j g consequently & 

residual anistropy factor (here cosh(K£>a)) between the 
two directions, which docs not disappear at large dis- 
tances r from the charges. 



A *(r) 



2q cosh(rei}a) exp(-ftj)r) 
4-ke r 



2q exp( — K D r) 
4tt£ r 
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FIG. 1. Illustration of the anisotropic effect. In the x 
and 2 directions, the far-field potentials differ from a factor 
cosh(Koa) which does not vanish at any distance. 

The same result is expected to hold for anisotropic 
macromolecules, with a residual and potentially strong 
anisotropy at large distances. The corresponding general- 
ization of the DLVO theory is thus required. We empha- 
size immediatly that the proposed description is mostly 
relevant in the case of moderately dissymetric objects, 
i.e. not too large aspect ratio, since the interaction en- 
ergy we shall compute is valid for distances between the 
objects larger than their typical size (this precludes infi- 
nite objects). This is anyway the case for many macro- 
molecules (Laponite clays, Bochmite rods, etc.). 

The purpose of the present paper is twofold 



and a description of the edge effects on the cylin- 
ders. An approximate analytical model is proposed 
yielding results in good agreement with numerical 
calculations. Note that we chose the finite cylinder 
geometry, not only for its relevance for polyelec- 
trolytes, but also because we expect edge effects to 
be particularly marked. This geometry is therefore 
a " benchmark" for the study of anisotropic electro- 
static interactions. 

As in the original calculation of Verwey and Over- 
beck [1], the macromolecules are specified by a constant 
electrostatic potential on their surfaces and the electro- 
static potential in the electrolyte solution is described at 
the level of the linearized mean-field Poisson-Boltzmann 
equation. However we will show extensively in a sub- 
sequent paper [23] that this assumption is justified for 
colloids bearing a large constant charge on their surfaces 
[18]. For small surface charges, the sketch of resolution 
presented thereafter can also be easily adapted. 

This paper is organized as follows: 

• in the first sections, we present the general method 
we have developped to construct the solution of the 
problem. 

• we subsequently deduce the general formula for 
the interaction between two anisotropic colloids at 
large distances. This yields a formal expression of 
the anisotropic factors discussed above. 

• we then apply this general method to the spe- 
cific case of finite cylinders. We first obtain the 
charge distributions on the cylinder, exhibiting the 
so-called edge effects. The influence of electrolyte 
concentration and finite-size effects are discussed. 

• An approximate analytical model is eventually pro- 
posed to describe these effects, yielding results in 
quantitative agreement with the numerical solu- 
tion. 



II. GENERAL CONSIDERATIONS AND 
DESCRIPTION OF THE PROBLEM 

A. Method of resolution : the auxiliary surface 
charge 



• we shall first describe in a general way the far 
field interaction between two anisotropic macro- 
molecules. This will lead to a generalized DLVO 
interaction between two non-spherical molecules, 
with a formal expression of the anisotropic inter- 
action factor. 

• we shall then apply these results to the case of fi- 
nite cylinders. A byproduct of this part of the 
work is the charge carried by the finite cylinder 



We consider a single charged macromolecule embeded 
in an infinite electrolyte solution. The solution is char- 
acterized by a Debye screening length, Id = 1/kd and 
as emphasized above, we assume that the electrostatic 
potential at the surface of the macromolecule, <J>o, is 
held constant [1] . The electrical double layer around the 
macromolecule is described at the level of the linearized 
Poisson-Boltzmann theory. This relies on a mean-field 
description of the micro-ion clouds, together with a small 
potential assumption. An extensive discussion of all these 
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assumptions can be found in [1,18]. We anticipate how- 
ever that the assumption of a constant potential at the 
macromolecule boundary naturally emerges as an effec- 
tive condition to describe correctly the far field obtained 
within the full non linear Poisson-Boltzmann theory, for 
colloids with a large bare charge, provided Koa is not too 
small [18,23]. 

In this context, outside the macromolecule, the electro- 
static potential obeys the linearized Poisson-Boltzmann 
(LPB) equation 



A$(f) = K 2 D ${r) 



(2) 



together with the boundary condition on the molecule 
surface E 



$(f) = $ 



(3) 



Note that we assume that the macromolecule interior is 
empty of charges, so that $(r) = <& for an y point r inside 
the macroion (this amounts to write A$(r) = 0). 




FIG. 2. Geometry of the problem. A macromolecule C, 
with a surface potential $ = $o, is immersed in an infinite 
electrolyte. The permittivity of the macromolecule is assumed 
to be much lower than that of the solvent(water), so that the 
electrostatic potential is assumed to be constant in the inte- 
rior of C. 

The surface charge density, a, is then obtained 
from the derivative of the electrostatic potential at the 
molecule surface : 



(4) 



where n is the (outer) unitary vector perpendicular to E 

and the notation -— I stands for n ■ V<E> 
\dnj s+ 

The standard Green function formalism is too cum- 
bersome to be applied in its simplest version to solve the 
previous equations, Eqs. (2) and (3). This is due to 
the existence of a non- vanishing excluded region for the 
micro- ions (inside the macromolecule), where the LPB 
equation, Eq. (2), does not apply. In other words, the rel- 
evant Green's function for the problem depends on par- 
ticulc shape and size, which seriously limit its practical 



interest. To circumvent this difficulty, we have therefore 
introduced an auxilary system, in which the LPB equa- 
tion applies everywhere in the volume. This is defined 
as : 



for f £ C 
for rY C 
for r £ E 



A$(r) = K 2 D <S>(r) 
A$(r) = K 2 D $(r) 
$(r) = $ 



(5) 



The corresponding surface charge on the molecule, cr, is 
defined here in terms of the solution $f u ii(r) of the pre- 
vious system of equations : 



cr(f) 



d$fuii(r) 
dn 



g$full(g) 

dn 



(6) 



Of course, the solution of Eq. (5), $fuii(^), reduces to 
the solution of Eq. (2), ^ompty if) , outside the macro- 
molecule. This matching originates in the unicity theo- 
rem for the operator — A + nj^ with Neumann or Dirichlct 
boundary conditions (see [20] for a similar result concern- 
ing the bare Laplace operator A). 

Now, the solution of Eq. (5), $f u n(f) can be defined in 
terms of the surface charge a : 



$fuii(r) 



a{f) G(r, f 1 ) dS' 



(7) 



where G(r,r f ) is the screened electrostatic Green func- 
tion, Gir,?) = exp(-K£,|f - f'\)/(4ire\r - ?\). The 
unknown auxiliary charge, a, is found by inverting the 
boundary condition on the macromolecule. This can be 
explicitly written as: for any point r on the molecule, 



$ = JJ a(?)G{ry)ds' 



(8) 



The overall result of these general considerations is a 
formal solution of the LPB equation, Eq. (2), for any 
point outside the macromolecule: 



$(r) = J J air 1 ) G(f, f 1 ) dS' 



(9) 



with the auxiliary charge a defined in Eq. (8). 

To get back to the "real" charge on the macromolecule, 
one has to compute the surface charge density as a func- 
tion of the auxiliary quantity, a. Using the definition 

on any point r on the colloid sur- 



cr(f) = - 

face E, one obtains: 



dn 



cr(f) 



d[G(f,r)] 



dn 



dS' 



(10) 



£ + 



In practice, the calculation of a which requires the in- 
version of the boundary condition, Eq. (8), can be per- 
formed analytically for simple geometries only, spheres or 
infinite rods (see below). For more complex case, such as 
finite cylinders as considered in this paper, a numerical 
calculation has to be performed to compute the inverse 
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matrix of G(f, r') on the (discretized) macroion. We shall 
show however that a simple model can be proposed which 
yields results in quantitative agreement with the numer- 
ical calculation. 

In the case of a given surface charge cr(r), this method 
of resolution can also be used to calculate the electro- 
static potential $(r) outside the colloid by computing 
the auxiliary charge at any point r 1 of £ using Eq. 
10 and then applying Eq. (9). Once again, the auxiliary 
charge a is the most relevant parameter to deal with the 
electrostatic potential created by a colloid immersed in a 
ionic fluid. 



B. The spherical case as an illustrative example 

Before going further, we come back to the simple spher- 
ical problem, where all previous different quantities, such 
as the bare and auxiliary surface charge, can be explicitly 
computed either by solving the LPB equation straightfor- 
wardly, or by using the auxiliary charge method sketched 
above. 

We consider an empty sphere £ of radius a, at a con- 
stant surface potentiel <I>o- On the one hand, the solution 
of the LPB equation is the usual Yukawa potential : 

p -K D (r-a) 

$(r) = $ a (11) 

r 

d$ 

The surface charge a is defined as a = — e— — (r = a) and 

dr 

is therefore given by: 

a = eK D fl + — J $o (12) 

On the other hand, the auxiliary problem described 
above consists in a sphere S filled with the elec- 
trolyte. Using the screened electrostatic Green function 
G(r, r*) = cxp(— Kb\t — r'|)/(47re|r — r'\), one may invert 
the integral equation, Eq. (8), to obtain the (uniform) 
auxiliary charge : 

a = eK D [l + coth(K D a)] $ (13) 

It is then straightforward to show that performing the 
integral in Eq. (10) allows to recover the surface charge 
density obtained above, Eq. (12). 

This simple example illustrates the difference between 
the bare and auxiliary problems which we have intro- 
duced in the previous section and two ways to calculate 
the real charge a as a function of <f>o. The first method 
could only be used because we knew the formal solution 
of LPB for a sphere at fixed potential but this is an excep- 
tion rather than the rule. On the contrary, the auxiliary 
charge method, even if it seems less sraightforward in this 
case, is a systematic way to compute the solution of LPB 
for given boundary conditions. 

We now turn to the calculation of the interaction en- 
ergy between two macromolecules. 



III. FAR-FIELD INTERACTION BETWEEN 
ANISOTROPIC HIGHLY CHARGED COLLOIDS 

Before focusing on a specific geometry, we first use the 
previous results to describe the interaction between two 
anisotropic charged macromolecules. 

We consider two colloids Cj (i = 1, 2) separated by 
a distance r much larger than the typical size D of the 
colloids. As we already noticed in the introduction, it 
is important to note that the restriction r > D makes 
sense for moderately dissymetric macromolecules only. 
We assume at this level that the charge profiles <Ji(r), 
and equivalently di(r), are known. The position of each 
colloid d is characterized by fixing a (somewhat arbi- 
trary) origin Oi for the molecule (this may coincide for 
example with the colloid center if it is symmetrical) . On 
the other hand, we assume that the orientation of the 
anisotropic colloid is described a unit vector Ui pointing 
into a direction Qi and an angle ipi corresponding to a 
rotation of C around Uj. We finally define the colloid- 
colloid direction using the unit vector u = O1O2/IO1O2I 
and introduce the bisector plane, II, of [O1O2] and O 
the intersection of II with (OiC^)- It will proove useful 
to introduce of system of coordinates {0,x,y, z}, with 
the x-axis corresponding to the axis (O1O2) (see Fig. 3). 
The distance between O and a point P is denoted as p. 
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FIG. 3. Calculation of the electrostatic interaction between 
two dissymetric macromolecules. An arbitrary center Oi, 
a unit vector Ui and a rotation angle ifi are defined for 
each molecule. We denote as r — IO1O2I the distance be- 
tween the two molecules, while the unit vector u is defined 
as u — O1O2/IO1O2I. We eventually introduce the bisector 
plane IT and the intersection point O between II and O1O2. 

We shall estimate the interaction force (acting on one 
macromolecule due to the other) by integrating the elec- 
trostatic stress tensor, T, defined as [20] : 

=> ( eE 2 \ => 

T= {P+— \ I —eE ® E (14) 

where / is the identity tensor, E the electrostatic field 
and P the hydrostatic pressure. The force acting on the 
macromolecule C 2 can be written accordingly as 
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T dS 



(15) 



Note that the integral runs over the bisector surface II, 
and not the colloid surface. This is a consequence of the 

fact that the divergence of the eletrostatic stress tensor 

=► 

T vanishes outside the macroions. 

We emphasize that the following calculations are con- 
ducted in the far field limit where the distance r is larger 
than the Debye length to = up 1 - This will allow us to 
expand the various quantities in powers of p/r. No spe- 
cific assumption is done however on the ratio between the 
typical size of the macromolecule, a, and Id- 

Hydrostatic equilibrium and (linearized) Poisson- 
Boltzmann equations, respectively — gr&dp+ p E = and 



A$ = k 2 d $, allow to write P = P^ 



Note that 



the linearization of the PB equation is fully justified in 
the present case since in the far field limit (r 3> n D l ) the 
dimcnsionless electrostatic potential e^/ksT is expected 
to be small. One therefore obtains 



T= 



co + 2 €Kd 



$2 _ Igl 

2 



(16) 



We denote E a the component of E in the direc- 
tion a, a=x, y, z. Then, for P e II : E y , E z — 
0{p/r)E x . Therefore, E 2 = E 2 [l + 0{p 2 /r 2 )} and 
E a Ep - E 2 S af} = E 2 0{p 2 /r 2 ).^ 

This allows to rewrite the force F 2 acting on the colloid 
2 as 



Fi^yjj [n 2 D ^ 2 (p)-E x 2 (p)^dS 



(17) 



Both the potential $ and the electric field E x in this equa- 
tion can be estimated from the solution for the potential 
created by a single colloid, as obtained in the previous 
paragraph, as we now show. First Eq. (9) can be written 



exp(-K D \r-r'\) dS , 



47re|r — r' 



(18) 



For distances r much larger than the typical size a of 
the macromolecule Cj, one might expand the previous 
equation for small r' to obtain the leading large r contri- 
bution: 

= exp(- KD r) rr exp( _ KDU; , . p )dS > n 9) 

with u r = f/r. We introduce at this point the total 
auxialiary charge Zi = Jf E d{f)dS' and the angular dis- 
tribution fi(P) defined as, 



MP) 



(F)exp(-K D u r ■ r')dS' (20) 



Using these definitions, one gets eventually the electro- 
static potential at point P as 



MP) = 



fi(P)e- K °< 
Airer 



(21) 



At the order O(pfr), is is straightforward to check that 
one might replace u r by u in the anisotropic factor /j of 
the previous equation : /j only depends on the angular 
coordinates (characterized by Ui and <pi). Note that the 
dependance on ipi disappears for axisymmetric colloids. 
From now on, we will only consider such objects so that 
may write /, = fi(ui) for simplification. The potential 
created by colloid i therefore reads 



MP) 



Zj fi{uj)e~ 



(22) 



In the r 3> k. d x limit, the corresponding electric field re- 
duces to Ei — ± hd ®i(P) u, with a plus (rcsp. minus) 
sign for i = 1 (rcsp. i = 2). The total electrostatic po- 
tential $ on the mediator plane II is written as the sum 
of the contributions due to each colloids, $ = $i + <f> 2 : 



$(P) - [Z x h(ui)+Z 2 f 2 {u 2 ) 



Aire r 



(23) 



Note that the superposition assumption for the potential 
is justified in the far field limit, where one may neglect 
mutual polarization effects. The same holds for the elec- 
tric field : E x = E\ + E 2 , leading to : 



E x = n D (Z x /i(tTi) - Z 2 / 2 (t5)) 



47rer 



(24) 



Introducing these expressions into Eq. (17) yields the 
following expression for the force F 2 : 

_ 2k 2 d Z x Z 2 fi(ui)f 2 (u 2 ) 



(47r) 2 e 



Jo 



2tt p dp ■ 



-2K D y/d 2 +p 2 
d 2 + P 2 



(25) 



In the far field region, r 3> Kjj 1 , it is legitimate to expand 
the integrand in powers of p/r and keep only the leading 
order : using e~ 2KD ^+^ = r (W/^+of/A- 4 )^ 
one may compute the integral to get 

ZiZ 2 h(u-[)f 2 (u 2 )e- K ° r _ 

F 2 = n D u (26) 

Airer 

which is always repulsive [24]. This force derives form 
the potential energy (again at leading order in kdt) : 



U 12 (r) = 



ZiZ 2 /i(«i)/ 2 («2)e-^ 
47re r 



(27) 



This expression for the interaction energy between 
the two macromolecules is one of the main results of 
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this paper. This generalizes the DLVO calculation for 
anisotropic molecules. Note that, in view of the various 
expansions performed, this expression is valid in the far- 
field limit, i.e. for interparticles distances r larger than 
both the Debye length and the typical size of the colloid 
a (say, to fix the ideas, r > 4£u,4a). 

As anticipated in the introduction, the interaction does 
not reduce at any distance to the isotropic DLVO result, 
obtained for spheres. The anisotropy of the interaction is 
described by the angular distribution fi(ui) and f 2(1^2) 
defined in Eq. (20). The latter is defined in terms of 
the (auxiliary) charge distribution on the macromolecules 
cr(r), or equivalently as a function of the bare surface 
charge a(r) using Eq. (10). 

We conclude this part by showing that the previous 
expression for the interaction energy indeed reduces to 
the standard DLVO expression for spheres (as it should). 
In this case, the bare and auxiliary surface charge on 
one sphere have been computed in the previous section, 
in Eqs. (12) and (13). On the other hand, the angular 
factor fi for each sphere i can be easily computed and 
reduces to f; = smhKDa . The latter is of course indepen- 
dent of any angular variable. Gathering these results, 
one retrieves the DLVO expression, Eq. (1) : 



A. Sketch of the numerical method 



We now consider a cylinder C with radius R and length 
L, at a contant potential $o- The resolution first starts 
with the computation of the auxiliary surface charge by 
inverting Eq. (8). This calculation involves the Green 
function G(f,r l ), G(f,r') = exp(— nu\r — r'|)/(47re|r — 
r'|), expressing the potential at point f created by a 
unit point charge in r. However due to the cylindrical 
symetry of the problem, one might reduce the dimension- 
ality of the problem by integrating the Green function on 
a ring (or small cylinder) whose center matches the axis 
of symetry of the cylinder, as illustrated on Fig. 4. This 
specific problem is considered below. Once the corre- 
sponding reduced Green function is known, the numerical 
task simplifies into a standard inversion problem. First, 
the cylinder C is decomposed into the superposition of 
small cylinders (on the lateral surface) or rings (on the 



head surfaces), denoted as Cj and IZk, with dimension I 
and surface charge density di (see Fig. 4)). 



U 



12 



1 + K£>a 



Zi Z 2 e 



Airer 



(28) 



A final note concerns the case of colloids with vanishing 
internal volumes. In the latter case, the bare and aux- 
iliary charge coincide, a — a, and our calculation leads 
back to the expression found in a different context by 
Trizac et al. [17]: 



U 



Z\ Zi e" 



12 



47re r 



- fi (ui) h{u 2 ) 



(29) 



IV. 



CHARGE DISTRIBUTION ON A FINITE 
ROD-LIKE POLYELECTROLYTE 



We now use the previous results to predict the far field 
interaction between two finite rod-like polyelectrolytes. 
In contrast to the spherical case, briefly considered in 
the previous section, the surface charge cannot be ob- 
tained analytically in this situation. Therefore, we shall 
first obtain numerically the surface charge on the cylin- 
der, by solving Eq. (8). We will then propose a simple 
analytical model yielding an approximate surface charge 
in good agreement with the "exact" numerical results. 

We emphasize at this point that the finite cylinder ge- 
ometry should be merely considered here as a generic 
situation where end effects are important. The present 
description could be easily extended to other related ge- 
ometries, like sphero-cylinders, ellipsoids, etc., though no 
fondamcntal difference is however expected. 




FIG. 4. The numerical calculations are performed by de- 
composing the cylinder C into small cylinders Cj of radius R 
and height £, and in rings IZk of radii ro and of width I. Each 
of these elementary surfaces carry a uniform surface charge 
density di. The numerical calculations were performed with 
I < 0.05£ D . 

Then Eq. (8) is discretized according to the equation : 



VjeE, $(»=;•) 



;•) = *o 



(30) 



where Gi (t\ , fj ) is the electrostatic potential created on 
the cylinder Cj or ring IZj by the cylinder Ci or ring IZi, 
carrying a unit surface charge density. 



G 



B. Reduced Green function 



C. Calculation of the surface charge 



As mentionned above, the previous inversion requires 
the knowledge of the potential created by an elementary 
ring or cylinder, which we now compute. To this end, we 
make use of the explicit expression of the electrostatic po- 
tential created by a disk of radius R at heigh z' carrying a 
uniform surface charge density (here equal to unity) and 
immersed in an electrolyte with Debye length In- This 
expression can be found in Rcf. [21] and reads : 



G disk (R,r,z) = 



R f°° .h{kR) Jo(fcr) 



2c: 



x cxp 



z-z'|) dk (31) 



with Jo and J\ the Bessel functions of order and 1. 
This is namely the potential created by a disk with ra- 
dius R at a point M, with cylindrical coordinates {r, z} 
(the origin being placed at the center of the disk). Note 
also that the dimension of Gdisk is given by R/e, since 
Gdisk is the potential created by a unit surface charge. 

Now the potential dG cy i(R 7 r, z) created, at a point M, 
by an infinitesimal cylinder with height dz', radius R and 
unit surface charge can be deduced directly as : 



dG cyl (R,r,z) = dz 
This leads to 



,dG d i sk (R, r,z) 
dR 



(32) 



dG cy \(R, r,z) = 



Rdz f°° k Jo(fcr) J (kR) 



2e 



\/k 2 + "d 



x cxp (-\/k 2 + k d 2 \z - z'\j dk (33) 

where the identity -^—[x J\(x)] — x J (x) has been used. 
dx 

As a result, the electrostatic potential created by a 
cylinder of radius R, height I and unit surface charge, 
with a center located in (0, z'), is given by 



R 

G cyi (R,z',r, z) = — 



z'+t./2 



dz" 



k Jo(fcr) J Q (kR) 



z'-t/2 Jo \Jk 2 + K 2 D 



x cxp (-^/k 2 + k d 2 \z - z"\J dk (34) 

Along the same lines, the potential G r i ng (ro, r, z) cre- 
ated by the ring of radius ro and of thickness £ can be 
expressed in terms of Gdisk (ro, r, z) according to the re- 
lation 

Gring(ro, r, z) = G d isk(r + £/2, r, z) - G d i s k(r - t/2, r, z) 

(35) 

where Gdisk(-R, r, z) is given above in Eq. (31). 

Note that in order to avoid numerical problems, the 
previous integrals must be reformulated specifically for 
the case z — z'. 



Inversion of the equation (30) yields the auxiliary sur- 
face charge a. The "real" surface charge, a, can be de- 
duced from a using Eq. (10). In a discretized form, this 
reads : 



W6 £*™> 



(36) 



where Gi takes either the cylinder or the ring form, ob- 
tained in Eqs. (34) and (35). This equation involves 
various derivatives of the Green function at the cylinder 

surface, namely : r+ , (^) r=R+ 

and(^). 

It will turn useful to write all the results in terms of 
dimcnsionless variables. All the lengths (such 
or L) are expressed in units of the radius of the cylinder 
R : eg L adlm = L/R. In the same way, the electro- 
static potential $ and surface charge densities a become 
respectively $ adim = e$/fc B T and cr adim = ^l B Raje 
where we recall that Ib is the Bjerrum length defined 
by Ib — e 2 /(An e ksT) (for water at room temperature, 
£b = 7 A). We also introduce dimensionless Green func- 
tions, as G adlm — eG/R (see previous remark on the di- 
mension of G). From now on, the index "adim" will be 
ommitcd to simplify notations. 



D. Numerical Results 



The previous equations are easily implemented numer- 
ically, provided the various expressions of the Green func- 
tions are written in terms of well-converging integrals as 
mentionned above. 

To fix ideas the potential on the macromolcculc is as- 
sumed to be Vb — 100 mVolts, so that <&o = 4 (see how- 
ever Refs. [18] and [23] for further justifications of this 
choice). 



1. Surface charge profiles 



We now present the results for the surface charges on 
the lateral and the head of the cylinder, that we shall de- 
note respectively as cr cy i(z) and ahead (z)- We first focus 
on the shape of the profiles. 

Typical results for these profiles are shown on Figs. 5 
and 6. 
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FIG. 5. Plot of the reduced surface charge on the lateral 
side of the cylinder, \n[a C y\(z/R)/a cy i(0) — 1]. The aspect ra- 
tio of the cylinder is L/R = 20 and the screening factor is 
kdR = 1.0. Note that z is in units of the cylinder radius R. 
The solid line is the result of the full numerical calculations, 
while the dashed line is the result of the "four parameter" 
model described in the appendix. The dotted line is the (re- 
duced) auxiliary surface charge a cy \{z/ R)/a cy \(Q). Note that 
the edge effect spans over a smaller distance for the auxiliary 
surface charge, compared to the "real" charge. See text for 
details. 




FIG. 6. Same as Fig. 5, but for the surface charge pro- 
file on the head of the cylinder ln[<7head(.z/-R)/0"head(O) — 1]. 
On this figure, we have also plotted the predicted scal- 
ing for the divergence in the absence of salt kdR = 0, 
cr hc ad(2:/i?)/a hcad (0) = (1 - (r/R) 2 )- 1/3 (open circles). 

Qualitatively, the main striking feature of these pro- 
files is the diverging surface charge close to the edges of 
the cylinder. This is of course the well-known edge ef- 
fect which is expected for charged objects with uniform 
potential. In the absence of electrolyte (kd = 0) , the di- 
vergence of the surface charge in the vicinity of an edge is 
a classical result [20]. For an infinite conducting dicdrc 
with an edge angle /3, the surface charge density a is 



found to diverge in the vicinity of the edge as p^/P^ 1 
where p is the distance to the edge [20]. In the present 
geometry, corresponding to (3 = 3ir/2, the surface charge 
is expected to diverge as p -1 / 3 . For a charged object 
embedded in an electrolyte, i.e. kd ^ 0, the situation 
is more complex. However the divergence is expected to 
remain, as can be understood from a simple argument. 
As mentionned in paragraph II B, the surface charge on 
a sphere with radius a and constant potential $o reads 

a = eK D ^1 + -j^j $0 (sec Eq. (12)). Now using this 

relationship for a non spherical object, one finds that the 
surface charge a diverge at the points where the radius 
of curvature a vanishes. 

Figs. 5 and 6 show that the auxiliary surface charge 
a also exhibits an edge effect. However the latter is 
more localized close to the edge, compared to the "real" 
surface charge a. As for a, the divergence of a can 
be understood using the results for the sphere, a — 
£Kd[1 + coth(K£>a)] 3>0j which indeed diverges as the ra- 
dius of curvature a vanishes. However, the transition 
from a small a region to a large a region is much more 
marked for the auxiliary surface charge than for the 
bare charge. Indeed from the previous expressions for 
a and a, one gets cr(a) = a(a = oo) + 0(1/ kdcl), while 
a(a) = a(a = oo) + C(exp[— Kdo]). The large a limit is 
therefore approached much more quickly for the auxiliary 
charge than for the bare charge, which is in agreement 
with the stronger localization of the divergence of the 
auxiliary charge close to the edge. 

We now report in more details on the variations of 
these density profiles when the size of the cylinder L and 
the screening length varied. Generally speaking 

the geometry of the problem is characterized by two di- 
mensionless quantities : the aspect ratio L/R and the 
amounts of screening kdR- Some general trends for the 
surface charge profiles emerge when these quantities are 
varied. First, the lateral surface profiles is found to sat- 
urate as the aspect ratio L/R goes to infinity. On the 
other hand, the head profile is found to be barely depen- 
dent on the aspect ratio. One expects in fact that the 
cylinder length L will only play a role when it is smaller 
or equal to the Debye length, k^ 1 , say KpL < a, with 
a of the order of a few units to fix ideas. Therefore for 
a given screening KdR, the profile is expected to satu- 
rate for aspect ratio larger than L/R ~ cx/(kdR)- This 
rule of thumb is confirmed when kdR is varied. In the 
present study, we have verified this assertion in the inter- 
val kdR € [0.1; 1] (data not shown). Typically one finds 
a ~ 5. Finally it is interesting to compare both profiles 
with the edge effect divergence predicted in the kd = 
limit, as argumented above. Only the charge profile on 
the head is found to be in semi-quantitative agreement 
with this scaling, as shown on Fig. 6. Note that in order 
to symetrize the predicted divergence, we compare the 
head profile with a head (r) / a head (0) = (1 - (r/i?) 2 )" 1 / 3 . 
On the other hand this prediction is found to fail for the 
cylinder surface charge. This is expected since in most of 
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the present calculations, the length L of the cylinder is 
larger than the Debye length, so that the nr> = profile 
is only a very crude approximation. On the other hand, 
the radius of the cylinder is always smaller than the de- 
bye length considered, and for the head, the kd = 
profile should be a fair but not so bad approximation for 
k d R < 1. 



2. Total lateral and head charge 



A more global quantity of interest is the total charge 
on the lateral surfaces and on each head of the cylin- 
der, respectively denoted as Z\ at e and Zhcade (e being 
the elementary charge). It prooves in fact useful to con- 
sider the average surface charges on the lateral surface 
CT iat CraSC = ^iat/(27ri2£) and on each head of the cylin- 
der, <J^^ SC — Zhead/(7T-R 2 ) (note that we plot below 
the reduced surface charge densities introduced above as 
a* = AnisRcr/e). These quantities are plotted respec- 
tively on Figs 7 and 8 as a function of the length of the 
cylinder L/R for various screenings kdR- In the limit of 
large aspect ratio, both charges saturate to finite values. 
Moreover, both charges are found to be increasing func- 
tions of the screening kdR- This is expected, as can be 
understood from the spherical test case, Eq. (12), as a 
benchmark. 




FIG. 8. Same as Fig. 7, but with the charge on the head 

aver a} 
head 



of the cylinder, ff] avera « e 



ANALYTICAL DESCRIPTION OF THE 
SURFACE CHARGE 



In this section, we propose a very simplified descrip- 
tion of the electrostatic problem, which has the virtue to 
provide analytic estimates of the surface charges. This 
estimate will proove usefull in fine to compute the in- 
teraction between two rod-like polyelectrolytes. A more 
detailed approach, including a description of the edge ef- 
fect, is proposed in appendix A. 



A. Uniform head and lateral surface charges 




FIG. 7. Dependance of the averaged surface charge on the 



cylinder, <T^ ras ° , as a function of the aspect ratio L/R. The 
solid line is the result of the full numerical calculation, while 
the dashed line corresponds to the four parameter model de- 
scribed in the appendix. From bottom to top, the screening 
factors are KdR = 0.2, kdR = 0.5 and kdR = 1.0. The 
L = oo asymptotic values are in agreement with the analytic 
result, Eq. (42). 



We consider a "zeroth order" approximation of the 
problem, consisting in a cylinder with uniform charges 
on the head and on the lateral sides. More specifically 
we assume a uniform auxiliary charge profile. We denote 
<7 cy i and (Thoad the auxiliary surface charge on the cylinder 
and on the head, and by a cy \ and ahead the corresponding 
"real" surface charges. 

At this level of approximation, equation (30) relating 
the auxiliary surface charge to the potential <&o reduces 
to a 2 x 2 problem : 

$0 = CTcyl G cy \{R, R, 0) + 2 ahead GheadOR, R, L/2) 

$0 = a cyl G cyl (R,0,L/2) 

+ ffhead [G hea , d (R, 0, 0) + G head (R, 0, L)] (37) 

The surface charges on the head and on the lateral side 
of the cylinder are then obtained using Eq. (10) as 



CTcyl — Ccyl 
+ Ohcad 



dG cyl (R,r,0) 
dr 

dG hoad (R, r, L/2) 
dr 



R+ 
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Chead = 0" cy i 



+ tThead 



dG cyl {R,0,z) 

8z 



L/2+ 



f dGh ca d(R,r, z) 
V dz 



(38) 



o+ 



In the previous equations, the derivatives of the Green 
functions are expressed in terms of integrals of Besscl 
functions (see Eqs. (31), (34)), which have to be com- 
puted numerically for any L and k. The systems in Eqs. 
(37-38) can be easily inverted to obtain the expressions 
of cr cy i and (Thcad as a function of the aspect ratio L/R 
and screening kdR- 

We do not report here the full expressions. Rather we 
consider the asymptotic L — > oo limit, in which the sur- 
face charges reach finite values. Note that this limit is 
reached for sizes L larger than a few Debye lengths. 

In the infinite L limit, the various Green function may 
be computed, yielding 

G cyl (R,R,0) = I (k d R) K (k d R) 
G disk {R,R, L/2)=0 
G cyl {R,0,L/2) = K (k d R)/2 

G disk (i?,0,0) = (l-e-^ R ) /2k d R (39) 

In the same way : 

(d/dr)G cyl (R,r = R+,0) = k d RI (k d R) K^kdR) 



(d/dz)G cyl (R,0,z = L/2+) = 



{d/dr)G disk (R,r = R+,L/2) = 
(d/dz)G disk (R,0,z = 0+) = I 



The previous result for a™f olm correspond to the semi- 
infinite cylinder limit [18]. One may also verify on Fig. 
7 that this result does indeed match the L — > oo limit of 
the averaged cylinder profile. Note that in contrast, one 
may verify that the uniform surface charge on the head 



uniform 
head 



(40) 



is only a fair approximation to the numerically 
computed averaged surface charge, even in the L — > oo. 
This is because for the screening considered (kdR < 1), 
the head always feels the egde of the cylinder. 



B. Towards a description of the edge effect 

A simple extension of the previous modclization can 
be proposed : adding a " ring" on the edge of the cylin- 
der should allow to capture the main features of the edge 
effect. This can be done in a straightforward way, but 
the details of the calculation arc somewhat cumbersome. 
We therefore report the details of this approach in the ap- 
prcndix A. This "four parameters" model gives results in 
good agreement with the numerical solutions. This can 
be seen on Figs. 7 and 8, where the results of this model 
are displayed (as dashed lines) against the full numerical 
results. 

However, the interactions between two polyelectrolytes 
do not involve the "real" charge, but the auxiliary charge. 
As we show below, the results of the much simpler " uni- 
form" approach described in the previous paragraph will 
proove sufficient to describe the interaction between two 
rods. 



VI. INTERACTION BETWEEN TWO ROD LIKE 
POLYELECTROLYTES 



Gathering results, we obtain after inversion of Eq. (37) 



-uniform 

a cyl 



-uniform = ZnpR^O 
head ^ _ e - KD R 



$0 

Iq(k d R) K {k d R) 



1 - 



2I (k d R) 



(41) 



We now denote these profiles as "uniform" to avoid any 
confusion with the numerical results. Using Eq. (38), 
one gets the "real" surface charge densities : 

uniform ^ K D RKi(n D R) 
^cyl = $ 77-, 

y Ko(kdR) 
e- KDR $ 
2I (k d R)K (k d R) 



uniform 

^hcad — 



+ 



\ _ e -K>L>R 



1 - 



2I {k d R) 



(42) 



In the limit of large upR, the surface charges are linear 
in KdR- This is expected since in this limit, one retrieves 
the planar results for which a oc k$ - 



We eventually turn to the description of the interac- 
tion between two rod-like polyelectrolytes. Our starting 
point is the potential energy obtained in section III, Eq. 
(27). The two crucial ingredients in this interaction en- 
ergy are : the total auxiliary charge Z on the cylinder; 
and the anisotropic term, f(P), defined in terms of the 
auxiliary charge profile in Eq. (20). We recall here this 
expression : 

f(P) = l/Z Jj a(f<)exp(-K D u r ■ r')dS' (43) 

These ingredients can be therefore easily computed from 
the full numerical solution, once the auxiliary surface 
charge has been computed. 



A. Total auxiliary charge 

We show on Fig. 9 the size dependence of the total 
auxiliary charge Zis/R, for various screenings ke>R. As 
can be seen on this figure, the charge is mainly linear in 
L. 
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FIG. 9. Total auxiliary charge Ztoxts/R as a function of 
the size of the cylinder L/R. The solid line corresponds to 
the full numerical resolution, while the crosses are the result 
of the uniform model. The dashed line is the result of the 
four parameter model detailed in the appendix. The dot- 
ted line corresponds to the uniform model with finite L (see 
text for details) From bottom to top the screening factors are 
KdR = 0.2, k d R = 0.5 and k d R = 1.0. 

This result is compared with the predictions of the sim- 
plified models we have proposed in the previous sec- 
tion. Within the simple uniform surface charge model 
described in Sec. VA, the total auxiliary charge reads 



Z = 2ir R L a. 



-uniform , rj p2 - uniform 
cyl "T £~ It a head 



(44) 



Eqs. (41) reports the expressions of the reduced auxil- 
iary charges (recall that in the previous section, reduced 
variables have been used a* = AkIbRo y 'e). This leads 
eventually to the following expression of the total aux- 
iliary charge function of the aspect ratio L/R and 
screening ke>R : 




1 



2 R I (k d R) K (k d R) 
kdR 



1 



-kdR 



i 



i 



2I (k d R) 



(45) 



This prediction is plotted as crosses on the previous fig- 
ures, showing a relatively good agreement with the "ex- 
act" numerical results. The agreement might be slightly 
improved by considering the complete L dependence, 
while staying within the uniform model. This corre- 
sponds to solving the 2x2 system of equations, Eqs. (38), 
with a numerical estimate of the Green functions for fi- 
nite L. We have plotted the results of this approach as 
dotted lines on Fig. 9. This improves slightly the agree- 
ment especially for small L and n^R. We also present 
the results obtained using the "four parameters" model, 
described in the appendix. This model adds to the result 
in Eq. (44) the contribution of the rings which capture 
the edge effects. This model is not analytic either and as 



can be seen on Fig. (9), it does not improve much the 
agreement. 

We conclude here that the very simple analytic ex- 
pression in Eq. (44) provides a useful and trustworthy 
approximation for the total auxiliary charge which enters 
the interaction energy, Eq. (27). 



B. Anisotropic Terms 



We report on Figs. 10 and 11 the numerical results for 
the anisotropic terms f(P) for two cylinder sizes L/R = 8 
and L/R = 20. These functions have been obtained af- 
ter numerical integration of Eq. (43) using the numerical 
result for the auxiliary surface charge. On these figures, 
the anisotropic terms are plotted as a function of the tilt 
angle 6, between the axis z of the cylinder considered and 
the unit vector u linking the two cylinder centers (see e.g. 
Fig. 3). 
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FIG. 10. Plot of the anisotropic factor of the finite cylin- 
der, /(#)/ < / >, as a function of the tilt angle. The solid 
line is the result of the integration of Eq. (43) over the numer- 
ically computed surface charge on the cylinder. The circles 
are the result of the uniform model (see text for details) while 
the dashed line is the result of the four parameter model de- 
scribed in the appendix. The aspect ratio is L/R — 8 and the 
screening factors are KdR = 0.2, 0.5 and 1.0 (from bottom to 
top for = 0). 
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FIG. 11. Same as Fig. 10 but for an aspect ratio L/R = 20 

It is instructive to compare these "exact" anisotropic 
factors to the predictions of the simplified models for the 
surface charges discussed in the previous section, Sec. V. 
Again, let us first concentrate on the uniform (auxiliary) 
charge model, proposed in section V A. In the frame of 
this simplified description, the anisotropic factor, in Eq. 
(43), can be computed analytically since the auxiliary 
charges are constant over the head and the lateral side 
of the cylinder. This leads the following expression for 

f(0) ■■ 



f(n) 



^/cylW + ^/headW 



(46) 



where Z cy \ — 2ir R L a cy \ is the total charge on the lateral 
sides of the polyelectrolyte, and Zhcad = 2tt R 2 <7 n0 ad is 
the charge on the heads of the polyelectrolyte. Using ex- 
pressions, Eqs. (41) obtained within the uniform model, 
one has 



J cyl R 



= - — $0 



1 



7, 

^ncad u 



2R I {k d R) K {k d R) 
kdR^o 



1 



2Iq{k d R) 



(47) 



and the total charge Z is given in Eq. (45). On the other 
hand, the expressions for anistropic factors due to the 
cylinder and due to the heads read 



/ cy i(0) = Iq(kdR sin 9) 



sinh 



' do L cos ( 



( hdL cos 9 



2I 1 ( KD RsmO) 

/hcad(6>) = p . a COsh 

n D R smti 



— ) 

kdL cos 8 



(48) 



This expression for /(#), using the previous expressions 
for Z cy \ and Zhcad, i s plotted against the numerical re- 
sults on Fig. 10 and 11 for two aspects ratios (L/R = 8 
and L/R = 20 respectively). The agreement is seen to 
be surprisingly good in view of the simplicity of the mod- 
elization. 



On these figures, we also show the prediction of the 
more detailed "four parameters" model, which includes 
a crude description of the edge effect, as detailed in ap- 
pendix A. This approach adds a contribution from the 



rings to the previous anisotropic factors, 



/ring(^) 



where Z, 



ring 



2 x 2irR£(a 3 + a 4 ) 



z J ring V 

is the total charge on 



the rings (see appendix A for details). The contribution 
to the anisotropic factor due to the ring, / r i ng , reads ex- 
plicitly : 



/ring(fl) = Io(kdR sin 6) cosh ( ^ DL ^° h ® 



(49) 



As can be seen on the figures, this more detailed descrip- 
tion does not improve much the agreement compared to 
the much simpler "uniform" approach. 

Such a good agreement using a very simple descrip- 
tion of the surface charge calls for some comments. The 
crucial point is that the interaction energy involves the 
auxiliary charges and not the bare charges. The full nu- 
merical resolution shows in fact that the edge effect is 
much more marked for the auxiliary charges than for the 
"bare" charge, in the sense that the divergence of the 
surface charge occurs much closer to the edge for the 
auxiliary charge. We have discussed this effect in section 
IV D 1. As a result, the auxiliary charge profile is more 
flat than the "real" charge profile. This feature allows to 
understand why the uniform model yields results in good 
agreement with the numerical results for the anisotropic 
factors. 



VII. CONCLUSION 

In the present paper, we have proposed a framework al- 
lowing to generalize the DLVO interaction for anistropic 
macromolecules. The central result is the electrostatic in- 
teraction energy between two anistropic macromolecules 



U 12 (r) = 



ZiZ 2 h{ui)f2{u2)e~^ r 

47T6 r 



(50) 



The main point resulting from Eq. (50) is that in a 
medium with finite salt concentration, the anisotropy is 
remanent at all distances. We have quantified this ef- 
fect and obtained general formulae for the anisotropic 
factor f(u, ip) (which only depends on u for axisymmet- 
rical objects) in Eq. (20). We have then applied this 
framework to finite rod-like cylinders. The previous cal- 
culations provide a simple and efficient description of the 
interaction between two such polyelectrolytes. In par- 
ticular, the simple uniform model leads to an analytic 
expression for the total auxiliary charge and anisotropic 
terms which enter the interaction energy, that turn out 
to be in good agreement with the full numerical solution. 
With this approximation, the anisotropic factor ,f{u) for 
a finite-size cylinder of length L and radius R at fixed 
potential $ takes a simple form 



12 



m = ^ff cyl (e) + ^f head (e). (5i) 

In the above expression, the auxiliary charges Z cy \, ^head 
and Z, as well as the anisotropy factors f cy i{@) and 
/head(0) are given by 



1 L 



1 



Zcyi'g 2 I ( KD R) K ( KD R) 



7, J^L- 
^ncad n 

z = 



Z C y\ + -^head 



1 



1 



2I (k d R) 



and 



sinh ( KdL cos0 ) 
f cy i(6)=I {K D Rsm6) ' ; 



( 



kdL cos 9 1 



/head (6*) 



2h{n D R sing) 
sing 



cosh 



kdL cos I 



As will be shown in [23] , the above expressions with the 
relevant choice of <&o almost corresponds to the interac- 
tion energy of two highly charged colloids far away from 
each other, irrespective of their bare charge. 
A few further comments are in order: 

• First, the interaction energy, at a fixed center to 
center distance between the two cylinders, is found 
to be minimum when the tilt angle (made between 
each cylinder and the center to center direction) is 
equal to n/2 i.e. when both cylinders axis are per- 
pendicular to the center to center vector. Apart 
from that, the angle between the two axis of the 
cylinders is not constrained, at this level of approx- 
imation (the two axis may equally be perpendic- 
ular or parallel). This is a consequence of retain- 
ing only the leading order contribution in the po- 
tential, and higher order terms (in exp(— K,£,r)/r l 
with i > 1) would split the aforementioned degen- 
eracy, and clearly stabilize the crossed rods com- 
pared to the parallel situation. On the other hand, 
the interaction is maximized when the two rods are 
coaxial (vanishing tilt angle). This result somehow 
contrasts with the infinite rod situation [14], for 
which the minimum energy situation corresponds 
to crossed rods (which is compatible with what we 
found), but with a totally different angular depen- 
dence, and also a different distance dependence. 

• The anisotropic term in the interaction potential 
results in a coupling between orientational and 
translational degrees of freedom. The strength of 
this anisotropy is moerover found to increase with 
salt concentration. These ingredients suggest that 
at high salinity, frustrated phases might form, in- 
dependently of van der Waals forces. However a 
full exploration of the phase diagram of charged 
rods using these previous results is required before 



reaching a definite conclusion on the formation of 
gels in rod like systems at large salt concentrations, 
as seen experimentally [6,8,9]. 

Work along these lines is in progress. 



Acknowledgments: We would like to thank Miguel 
Aubouy for inspiring discussions and an enjoyable col- 
laboration on related topics. 



APPENDIX A: A SIMPLE DESCRIPTION OF 
THE EDGE EFFECT 

1. general framework 

In this appendix a more detailed description of the 
edge effect is proposed. We extend the model described 
in section VA by incorporating a specific charge on the 
edge of the rod-like macromolecule. More specifically, we 
model the auxiliary surface charge as the superposition 
of a uniform charge on the head and on the lateral sur- 
face of the cylinder, supplemented by a ring charge on 
the edge of the macromolecule, as shown on Fig. 12. 




FIG. 12. Simplified description of the edge effects 

From a technical point of view, we separate the ring 
charge on the edge of the molecule as a ring of radius R 
on the head, and a ring of radius R on the lateral side of 
the cylinder (see figure 12). The extension of the lateral 
ring is denoted as £ cy \, and that on the edge £head- There 
is therefore four parameters in the model: respectively 
the uniform surface charge on the head dhcad, on the lat- 
eral sides er cy i, <T cy i. d gc and <7h ea d edge ■ In the following 
results, we have chosen i cy \ = £hcad = 0-05R. Results are 
only weakly dependent on this choice. As in section V A, 
one has to solve Eq. (30), relating the auxiliary surface 
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charge to the potential $o- Within the simplified analy- 
sis, and taking into account the symetry of the cylinder, 
this equation reduces to a 4 x 4 inversion : irrespective 
of j 



(Al) 



where the summation i runs overs the different parts of 
the simplified object : eg, i = 1 stands for the center of 
the heads of the cylinder; i = 2 stands for the middle part 
of the cylinder; while i = 3 and i = 4 stand for the rings 
on the edges. The "Green functions" Gij are defined in 
terms of the Green functions Gdisk, G cy i and G r mg whose 
expressions are given respectively in Eqs. (31), (34) and 
(35) (see below for the detailed expressions of the 4x4 
matrix Gij). Once the auxiliary charges &i are known, 
one obtains the "bare" charge a everywhere on the cylin- 
der using Eq. (10). This can be written formally : 



. dG{r, i) 



dri 



(A2) 



where the notation G(r, i) stands for the Green function 
computed at point r due to charge i defined above; d/dn 
denotes the derivative along the normal to the surface at 
point r. 

The previous equation, Eq. (Al), is easily inverted and 
the corresponding surface charges are plotted on the pre- 
vious figures, Figs 7 and 8. As shown on these plots, the 
approximate description yields results in excellent agree- 
ment with the full numerical resolution for any aspect 
ratio L/R and screening kdR- 

As a consequence, despite its simplicity, the simplified 
description of the auxiliary charges contains most of the 
physics of the edge effect. Also, as shown by the previous 
argument, a better agreement is expected for large KrjR. 



2. technical details 

The cylinder is decomposed into 4 different pieces: 

• a lateral part of length L and radius R 

• two disks of radius R 

• two lateral rings of radius R, of hcigths 4 y i and of 
centers located in [0,±(L/2 =F 4yi/2)] 

• two rings of radius R — 4cad/2 and widths 4ead 

respectively denoted 1,2,3 and 4 and carrying the 
uniform surface charge densities a cy \, <7h oa d, ^cyi edge 
and ahead edge • To simplify the formulation of the 
equations, we respectively call G cy i(r , r, £, z) and 
Gring^Oj 4cad, f, z) the electrostatic potentials by a cylin- 
der of radius ro and of height £ located in (r, z) and by 
a ring of radius r and of width 4ead in (r, z) with the 



origin of the coordinates (0, 0) located in the center of 
the cylinder or of the ring. 

In order to find the auxiliary charges on the disks, rings 
and lateral sides of the cylinder, one has to solve the 4x4 
linear problem, obtained from Eq. (30) : 

Vje {1,2,3,4},^^^- =$ 



The coefficients Aij are given in terms of the expres- 
sions of G cy i and Gdisk given in Eqs. (34) and (35) : 



A n = G cyl (R,R,L,0) 

A 12 =2G disk (R,R,L/2) 

Ai 3 = 2 G cyl (R, R, £ cyh L/2 - £ cyl /2) 



(A3) 



Ai4~ 2 G Ilng (R — €head/2, R, L/2) 

A 21 =G cyl (R,0,L,L/2) 

A 22 = Gdisk (R, R, 0) + Gdisk (R, 0, L) 

A 2 3=G cy l(i?,0,4yl,4yl/2) 

+G cyl (R,0,£ cyh L-£ cyl /2) 

A 2 4=G ri „g( J R-4ead/2,0,0) 
+ G r i ng (i? — 4ead/2, 0, L) 



(A4) 



(A5) 



A 3 i=G cyl (.R,i?,L,L/2-4yi/2) 

A32= Gdisk(-R, R, 4yl/2) + Gdisk(-R, R, L — 4yl/2) 

A 33 = G cyl (R, R, 4 y i/2, 0) + G cyl (R, R,L~ 4 y i) (A6) 

^34= G r ing(i? — 4oad/2, R, 4yl/2) 

+G r ing(-R — 4oad/2, R,L — 4yl/2)) 

A 4 i= G cy \(R, R — 4oad/2, L, L/2) 

Ai2= Gdisk(-R, R — 4cad/2, 0) + Gdisk(-R, R — 4ead/2, L) 

A 4 z= G cy \(R, R — 4 y i/2, 4yi/2) + G cy i(R, R — 4 y i/2, L — £ cy \/2) 

A 44 = Gdisk(-R - 4oad/2, R — 4oad/2, 0) 

+G d isk(i? - 4cad/2, R - 4cad/2, L) (A7) 

Once the a have been calculated, we get cr using Eq. 
(36), which reads within the simplified description : 



V* € {1, 2, 3, 4}, Gi{r) = J2 *i Bij (r) 

3 

The coefficients (r) are given by 
'dG cy i(R, r,L,zY 



(A8) 



Bn(zy- 



B 12 (z)= 



+ 



dr 



f dG disk (R,r,L/2-z) 
\ dr 
gGdiskCg, r, L/2 + z) 
dr 



R+ 
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+ 



B 14 (z)-- 



+ 





dr 


dG ,(R r 


£ i 1,11 — P 1/9 + 7"! 

<-cyl )-'-'/ ^ <-cyl/ z T 6) 




dr 


dG r i ng (R - 


-4ead/2,r, L/2-z)\ 




dr ) i 


dG I i ng (R - 


-£ head /2,r,L/2 + z)\ 



(A9) 



R+ 



R+ 



B 2 i(r)-- 



dr 

dG cyi (R, r,L,z) 
dz 



R+ 



R+ 



L/2+ 



B22{r)= l + fdG disk (R,r,z) 



dz 



L+ 



B 2 3{r)-- 



+ 



B 2i {r) = 



dG cyl (R,r,£ cyl ,£ cyl /2) 
dz 

dG cy i(R, r, £ cyh z) 
dz 



(A10) 



(L-£ cyl /2) + 
Iring(r) , ( dG, ing (R - 4cad/2, r, z) 

dz 



(L-£ cyl /2)H 



with l r ing(7) = 1 if R — 4cad < t < R and otherwise. 
fdG cyl (R,r,L,L- 4 y i/2) s 

^32 = 



9r 

dG dlsk (R, r,£ cyi /2) 



R+ 



dr 



R+ 



+ 



B33 — 



+ 



dG disk (R,r,L-£ cyl /2) \ 

dr )r+ 
( dG cyl (R,r,£ cyl /2,0) \ 

V dr J R+ 

dG cy i{R, r, 4 y i, L - £ cy \) 



(All) 



dr 



R+ 



B- 



34 = 



R+ 



f dG ring (R-£ hcad /2,r,£ cyl /2) 
; . dG ring (R - 4cad/2, r,L- £ cyl /2) \ 

dr ) R+ 



Ba-i — 



dG cyl {R,R-£ head /2,L,z) 
dz 



L/2+ 



1 /<9Gdisk(^, -R — 4ead/2, L) 

Bi2= 2 + I 9i 



43 = 



/ ggcylCg, ~ 4cad/2, 4yl, 4yl/2) 

I 9z 



(A12) 



<cyl/2-» 



+ / 9Gcyl(^, - 4cad/2, 4yl, L - 4yl/2) \ 

dz ) (L-e cyl /2)+ 



Baa — 



1 / dG T i ng (R — 4cad/2, i? — 4ead/2, i 



<9z 



L + 
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